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1.  Introduction 


In  conventional  vehicles  using  internal  combustion  engines,  approximately  40%  of 
the  fuel  chemical  energy  available  is  lost  to  the  exhaust  stream.  Much  effort  has 
been  directed  at  recovering  some  of  this  waste  heat  for  conversion  to  useful  energy.^ 
One  of  the  most  commonly  proposed  methods  for  vehicle  waste  heat  recovery 
(WHR)  is  the  use  of  a  thermoelectric  generator  (TEG),^’^  although  other  techniques 
have  been  suggested  including  thermodynamic  power  cycles'^  and  direct  thermal 
coupling.^  Thermoelectric  (TE)  materials  offer  several  advantages  that  make  them 
particularly  attractive  for  mobile  applications  including  being  lightweight,  solid- 
state,  and  passive.  This  has  the  potential  for  producing  a  reliable,  low  complexity 
solution  for  the  vehicle;  however,  current  TEGs  suffer  from  low  conversion 
efficiencies,  especially  at  elevated  temperatures  (exhaust  gases  can  range  from 
800-1100  K).^  This  creates  a  significant  engineering  challenge  in  integrating  the 
TE  to  the  exhaust  system  and  identifying  a  useful  configuration  with  respect  to 
recovered  power. ^  While  much  effort  is  going  into  the  improvement  of  TE  material 
performance  and  several  studies  have  shown  the  potential  of  TE  materials  with  a 
dimensionless  figure  of  merit  ZT  >  1,  delivering  these  materials  outside  the 
laboratory  environment  has  proven  challenging.^’^  As  such,  implementing  a  TEG 
for  vehicular  exhaust  energy  recovery  requires  careful  attention  be  paid  to  the 
thermal  integration  of  the  TEG  so  as  to  take  maximum  advantage  of  the  conversion 
efficiency  of  available  materials. 

When  designing  a  TE  WHR  system,  care  must  be  taken  in  making  simplifying 
assumptions  regarding  both  vehicle  thermal  and  energy  conditions  and  TEG 
operation.  Because  of  the  inherently  low  conversion  efficiency  of  available  TE 
materials,  improper  assumptions  significantly  alter  the  perceived  benefit  of  a 
particular  integration  strategy.  Many  studies  evaluating  the  viability  of  TEG  WHR 
have  made  assumptions  regarding  the  operating  point  of  the  TEGs,  namely,  fixed 
temperature  conditions  at  the  TEG  hot  and  cold  sides,  constant  TE  properties,  no 
consideration  of  heat  rejection  to  the  coolant  loop  or  energy  cost  for  rejecting  that 
heat,  no  impact  of  exhaust  heat  exchange  on  engine  performance,  and  no 
accounting  for  circuit  losses  in  converting  the  electricity  to  usable  form.'*’"'^  In 
reality,  vehicle  exhaust  conditions  (temperature  and  enthalpy)  vary  widely  with 
engine  speed  and  load.'^  As  a  result,  recovery  systems  designed  around  the  TEG 
best-case  operating  point  achieve  efficiency  metrics  on  paper  but  fail  to  do  so  in 
practice,  because  real  applications  will  spend  little  time  at  any  particular  operating 
point,  much  less  the  optimal  one. 

One  possible  method  of  addressing  the  first  concern — spatial  temperature  non¬ 
uniformity  along  the  exhaust  stream — is  to  use  a  high  thermal  conductivity  device 
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(e.g.,  heat  pipes/vapor  chambers  or  heat  spreaders  with  integrated  diamond  layers 
or  carbon  nano  tubes  arrays)  to  load  level  the  temperature  of  the  TEG  hot  sides. 

The  advantage  of  this  would  be  more  uniform  power  generation  across  the  array  of 
TEGs  when  compared  to  a  configuration  where  TEGs  along  the  length  of  the 
exhaust  are  exposed  to  progressively  lower  temperatures,  which  shifts  the  material 
properties  and  leads  to  non-uniform  power  generation.  Additionally,  it  may  prove 
to  be  a  simpler  and  less  expensive  solution  than  trying  to  optimize  the  TEG  material 
configuration  to  match  the  expected  spatial  temperature  gradient  down  the  length 
of  the  exhaust. 

In  this  report,  a  heat  transfer  model  of  the  TEG  is  constructed,  and  the  impact  of 
using  high  effective  thermal  conductivity  devices  to  spatially  load  level  the  TE 
array  is  considered.  We  examine  the  effect  of  the  exhaust  stream  to  TEG  hot  side 
thermal  resistance  on  temperature  load  leveling  and  demonstrate  that  TEG 
temperature  profiles  approach  the  commonly  assumed  constant  temperature 
difference  condition  when  heat  spreading  devices  are  used. 

2.  Model  Formulation 


One  of  the  more  complete  system-level  TEG  WHR  models  in  literature  is  that  of 
Goncalves  et  al. not  defined,  order  to  facilitate  comparison  with  their 
work,  we  use  a  similar  TEG  configuration.  This  study  considers  one  16x6  array 
of  TEGs  (9505/127/150B,  Eerrotec,  as  shown  in  Eig.  1)  integrated  in  cross-flow 
heat  exchangers  operating  at  steady-state  and  subject  to  various  boundary 
conditions  (Table  1).  Traditional  (Eig.  2a)  and  heat  spreader  assisted  (Eig.  2b) 
arrangements  are  evaluated.  In  both  cases,  the  heat  flow  path  from  the  exhaust 
stream  to  the  TEG  hot  sides  is  described  by  a  thermal  resistance  Rhot,  which  is  varied 
to  capture  the  effect  of  a  range  of  heat  transfer  techniques  (e.g.,  surface 
enhancements,  fins,  interfacial  thermal  greases).  Eor  the  traditional  cross-flow  case, 
the  lowest  Rhot  evaluated  for  each  set  of  boundary  conditions  yields  7)/  ~  500  K, 
which  is  the  upper  limit  for  the  TEG  property  data  used.  Eor  the  heat  spreader 
assisted  cases,  the  lowest  Rhot  evaluated  (0.01  K/W)  corresponds  to  an  idealized 
“best  case”  scenario,  as  achieving  it  would  require  multiple  densely  firmed  heat 
pipes,  a  high  forced  convection  coefficient  (>250  W/m^K,  the  high  end  of  the  range 
suggested  by  Incropera  and  DeWitt^°),  and  no  interfacial  and  spreading  resistances. 
On  the  coolant  side,  a  constant  thermal  resistance  from  the  coolant  stream  to  each 
module  is  taken  as  Rcooi  =  0.625  K/W,  which  corresponds  to  a  convection 
coefficient 

h  =  1000  W/m^K  and  a  heat  transfer  area  ^  =  1.6  x  10“^  m^  equal  to  the  footprint 
of  a  single  module. 
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Fig.  1  TE  generator  module  (9505/127/150B,  Ferrotec)  properties  as  a  function  of 
temperature.  The  lines  are  a  second-order  polynomial  fit  to  data  (from  Goncalves  et  al.^"' 
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Table  1  Exhaust  and  cooling  loop  conditions 


mexh  (kg/s) 

Te.h.i(°C) 

nicooi  (kg/s) 

Tcool,o  i°C) 

case  1 

0.02 

650 

0.15 

65 

case  2 

0.01 

580 

0.15 

65 

case  3 

0.01 

580 

0.15 

50 

case  4 

0.02 

650 

0.15 

50 
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CV 
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Fig.  2  a)  Traditional  and  b)  heat  pipe  assisted  cross -floss-  configurations,  c)  Control 
volumes  (CVs)  used  in  the  analysis. 

^^exh ^cool 


The  heat  spreader  assisted  aiiangement  employs  an  integrated  heat  pipe  and 
spreader  assembly  for  temperatiue  load  leveling  on  the  TEG  hot  sides.  Due  to  the 
high  effective  thermal  conductivity  of  the  heat  pipe,^^  a  hrmped  element  approach 
is  taken  (Bi  <0.1  for  heat  pipe  length  <2  m  assrrmiug  an  effective  thermal 
condirctivity  of  5000  W/mK),  and  the  heat  pipe  temperatiue  is  assumed  to  be  a 
constant  value  detennined  by  a  separate  iterative  calculation  that  matches  the 
exhaust  heat  extracted  by  the  heat  pipe  to  the  total  heat  conducted  tluough  the  TEG 
material  in  the  module  aiiay.  For  simplicity,  the  thennal  resistance  Rhot  accoimts 
for  all  resistances  between  the  exlraust  stream  and  TE  hot  sides  of  the  modules  in 
the  array.  Heat  spreaders  can  offer  effective  thermal  conductivities  of 
10000  W/mK;^^  thus,  all  TEG  hot  sides  are  assumed  to  be  mamtained  at  a  constant 
temperatiu  e  by  a  heat  spreader. 
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Energy  balances  capturing  the  Seebeck  effect  and  conduction  through  the  TEG 
modules  are  performed  on  the  control  volumes  illustrated  in  Eig.  2c.  A  system  of 
equations  is  formulated  (Eqs.  1-4)  and  solved  using  a  Newton-Raphson  iterative 
scheme. TEG  hot  and  cold  side  temperatures  are  assumed  to  be  constant  for  each 


module,  and  all  TEG  module  properties  are  evaluated  using  the  average 
temperature  at  the  current  iteration. 

TEG 

^  ( 

Rfiot 

'Texh,i+Texh,o  -  ^  S  1  +  K(Tff  -  Tc)  -  I^RjeI^ 

(1) 

1 

Rcool 

{Tc  -  =  S I  Tc  +  K{Ti,  -  Tc)  +  I^Rte/2 

(2) 

1  ( 

'P exh.i'^'P exh,o  )  T^exhCp.exhiT exh.^P exh,o) 

(3) 

Rfiot 

V  2  - 

1 

( rji  'Rcool,i'^'Rcool.o\  ^cool^p.coolir  cool, o~'R  cool, i) 

(4) 

Rcool 

2  )- 

S  and  Rte  are  the  module  Seebeck  coefficient  and  electrical  resistance,  r  is  the 
number  of  rows  in  the  TE  array,  and  Cp  is  taken  to  be  1000  J/kg  for  the  exhaust  gas. 
The  coolant  Cp  is  evaluated  at  each  iteration  based  on  a  50:50  water: ethylene  glycol 
mixture;^"^  the  coolant  mass  flow  rate  rh^ooi  is  taken  as  0.15  kg/s,  equal  to  the 
midpoint  of  flow  rates  evaluated  by  Oliet  et  al.^^  The  TEGs  in  the  array  are  all 
connected  in  series;  thus,  the  generated  current  /  must  be  equal  in  all  modules.  This 
adds  an  additional  constraint  to  the  system  of  equations: 


/  = 


(Th-  Tc) 


^load 

n 


+Rte 


(5) 


where  n  is  the  total  number  of  TE  modules  and  Rioad  is  load  resistance,  assumed  to 
bs  Rioad  =  El  Rte- 

The  power  generated  by  each  TE  module  is^' 

Pgen=Sl(.TH-  Tc)-1^/Rte  (6) 


The  TE  conversion  efficiency  rjconv  of  each  TE  module  is  evaluated  by  comparing 
the  power  generated  to  the  input  thermal  power:^^ 


Vconv  ~ 


Rgen 

K(Th-Tc) 


(V) 


The  system  conversion  efficiency  is  calculated  by  comparing  the  total  generated 
power  Pgen  to  the  exhaust  enthalpy  as 


r]  = 


liP, 


gen 


^exh^ppxhiTexh.i  P amb) 


(8) 


5 


where  Tamb  is  taken  to  be  300  K.  Also  of  interest  is  the  total  heat  for  rejeetion  on 
the  TEG  array  coolant  loop; 

Qc  ~  ^cool  ^P,cooi{TcooI,o  '^cool.i  )  (9) 

where  Tcooi.o  and  Tcootj  are  evaluated  at  the  array  outlet  and  inlet,  respectively. 

3.  Results  and  Discussion 


3.1  TEG  Analysis 

Figure  3  presents  the  hot  and  cold  side  temperatures  of  TEG  modules  along  the 
axial  length  of  an  exhaust  for  traditional  and  heat  spreader  assisted  cross-flow 
configurations  under  case  1  boundary  conditions.  Without  the  integrated  heat  pipe 
and  spreader  assembly,  the  TEG  module  temperatures  rapidly  decline  along  the 
length  of  the  exhaust  pipe.  The  temperature  non-uniformity  is  more  pronounced  at 
low  Rhot  because  more  heat  is  transferred  from  the  exhaust  stream,  leading  to  higher 
module  temperatures  as  the  coolant  outlet  temperature  remains  constant.  For  both 
configurations,  increases  in  Rhot  diminish  the  TEG  hot  and  cold  side  temperature 
difference,  which  adversely  affects  power  generation  (Eq.  6).  Dashed  lines  in  the 
figure  are  the  average  TEG  hot  and  cold  side  temperatures,  which  are  representative 
of  the  idealized  case  where  a  constant  AT=  Th-Tc  may  be  taken.  It  is  seen  that  the 
heat  pipe  assisted  configuration  provides  similar  temperature  profiles — Th  is 
constant  due  to  the  heat  spreader  and  the  rise  in  Tc  along  the  TEG  array  is  minimal 
due  to  the  high  coolant  flow  rate.  This  indicates  that  a  heat  pipe  assisted  cross-flow 
configuration  can  achieve  results  very  close  to  the  ideal  constant  zlT  case. 
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Fig.  3  T£  generator  module  hot  and  cold  side  temperatures  along  the  array  length  mth 
case  1  exhaust  and  cooling  loop  conditions.  The  dashed  lines  indicate  the  average  hot  and  cold 
side  temperatures  in  the  cross-flow  arrangement. 

The  effect  of  uou-imiform  temperatiue  on  the  TEG  properties  is  captiued  by 
plotting  ZT  along  the  length  of  the  exhaust  pipe,  as  shown  in  Fig.  4a.  Module  ZT 
decreases  sharply  with  increasing  temperatiue  (see  Fig.  1);  thus,  cases  with  higher 
Rhot  lead  to  higher  ZT  values.  The  generated  power  for  each  module  is  conelated 
to  S^{Th  -  Tcf  (Eqs.  5  and  6).  For  traditional  cross-flow  configmations,  the  first 
modules  in  the  array  are  subjected  to  the  highest  temperatiues,  and  AT  decreases 
along  the  length  of  the  array.  The  first  module  at  the  lowest  Rhot  (2  KAV  for  case  1) 
yields  low  power  generation  and  efficiency  due  to  its  very  low  value  of  S.  As  the 
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average  module  temperature  decreases  along  the  length  of  the  TEG  array,  the  power 
peaks  and  proceeds  to  decline  toward  the  end  of  the  array,  yielding  very  non- 
uniform  Unlike  the  traditional  cross-flow  configuration,  the  heat  pipe  assisted 
system  exhibits  a  zlT  that  increases  as  the  end  of  the  TEG  array  is  approached.  This 
is  due  to  the  fact  that  This  leveled  by  the  heat  spreader,  and  the  coolant  flow  enters 
the  heat  exchanger  at  the  end  of  the  array,  thus  yielding  the  largest  zIT  in  the  system 
(Eig.  4b).  TEG  conversion  efficiencies  are  below  1%  for  each  module  (Eig.  4c). 
Overall  device  efficiencies  are  even  lower  (Eig.  5a)  indicating  that  the  current 
system  does  a  poor  job  of  converting  available  exhaust  enthalpy  to  usable  power. 
Overall  efficiency  could  be  improved  by  adding  more  TEG  arrays;  however, 
extracting  more  heat  from  the  exhaust  stream  would  drop  the  TEG  hot  side 
temperatures  (and  thus  shift  the  operating  point  of  the  TEG  modules),  meaning  that 
the  results  of  the  present  investigation  cannot  be  directly  extended  to  a  multi-array 
system. 
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(a)  (b) 


(c) 


Fig.  4  TE  geaerxtor  module  a)  ZT,  b)  power  generation,  and  c)  conversion  efOriency  along 
the  array  length  with  case  1  exhaust  and  cooling  loop  conditions.  For  the  cross  flow  and  heat 
pipe  arrangements,  Ri,ot=  2,  5, 10, 50  KA\'  and  Rhot  =  0.01,  0.05,  0.1,  0.5  K/W,  respectively. 


3.2  WHR  System  Analysis 

It  is  worth  e.xamming  the  impact  of  the  TEG  WHR  system  on  the  rest  of  the  system 
to  better  imderstand  the  net  power  generation  provided.  First,  heat  thiough  the  TEG 
not  converted  to  useful  power  must  be  rejected  by  the  coolant  loop.  Figure  5b-c 
presents  the  total  system  power  generation  and  coolant  loop  heat  rejection.  Due  to 
the  low  TEG  conversion  efficiencies,  the  maximmu  ratio  of  Pgen  to  Qc  achieved  for 
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(%)  ^ 


both  cross-flow  and  heat  pipe  assisted  configiuatioiis  is  0.007,  iudicatmg  that  under 
the  conditions  of  case  1  (the  most  optimistic  scenario  considered),  1000  W  must  be 
dissipated  for  every  7  W  generated.  Taking  the  generic  radiator  studied 
parametrically  by  Oliet  et  al.^^  and  assuming  adequate  air  flow  over  the  radiator 
fins,  the  dissipation  of  10  kW  requires  about  3  W  of  coolant  pumping  power,  equal 
to  ~5%  of  Pgen  for  the  best-peiTonning  cases.  Tliis  does  not  accoimt  for  any 
additional  fan  power  that  may  be  required  to  generate  the  needed  auflow,  which 
may  be  significant. 


(a) 


(b) 


(c) 


Fig.  5  a)  Overall  efficiency  and  b)  power  generation  of  the  TE  generator  array  for  exhaust 
and  cooling  loop  conditions  1—4  as  a  function  of  Rhoi.  c)  Total  heat  for  rejection  on  the  coolant 
side  of  the  heat  exchanger. 
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A  secondary  system  concern  is  the  effect  of  the  hot  side  heat  exchanger  on  engine 
performance.  Previous  studies  have  shown  that  for  a  wide  range  of  engines 
additional  backpressure  can  decrease  engine  power  by  ~1%  per  inch  Hg.^^  A 
specific  exhaust  heat  exchanger  design  would  need  to  take  this  effect  into  account 
to  present  a  net  system  efficiency  increase.  Poor  heat  exchanger  design  has  the 
potential  to  negate  any  net  power  recovery  through  decrease  prime  power 
conversion  efficiency. 

3.3  WHR  Model  Comparison 

Goncalves  et  al.^®  built  a  computational  model  of  a  cross-flow  WHR  system  with 
the  TEG  hot  sides  mounted  to  a  solid  copper  heat  spreading  block.  Despite  the  fact 
that  modeling  results  showed  large  thermal  gradients  (87.5  K/m)  along  the  axial 
length  of  the  copper  block,  their  analysis  assumed  constant  hot  and  cold  side 
temperatures,  and  TE  properties  were  calculated  based  on  the  average  of  the  hot 
and  cold  side  temperatures.  While  calculating  the  appropriate  coolant  mass  flow 
rate  necessary  to  reject  the  estimated  39  kW,  the  authors  neglected  any 
consideration  of  thermal  resistances  in  the  heat  path  from  the  TEG  cold  side  to  the 
fluid  in  the  coolant  loop.  By  assuming  the  TEG  cold  sides  were  maintained  at  the 
average  coolant  temperature  (325  K),  they  essentially  assume  Rcooi  =  0  K/W,  which 
is  most  certainly  not  the  case.  Applying  these  conditions  and  assumptions  to  our 
model,  agreement  is  within  7%. 

As  noted  previously,  we  arrive  at  Rcooi  =  0.625  KTW  using  a  generous  convection 
coefficient  of  1000  W/m^K.  Eurther  reduction  in  Rcooi  is  non-trivial  and  will  add  to 
the  complexity  and  cost  of  the  coolant  loop.  As  a  result  of  this  assumption  in 
Goncalves  et  al.,'^  the  power  generation  and  efficiency  of  the  TEG  array  is  vastly 
overstated — they  estimate  the  power  generation  by  the  TEG  array  to  be  1530  W  at 
a  TEG  conversion  efficiency  of  3.9%.  Running  the  cross-flow  model  developed  in 
this  work  with  boundary  conditions  identical  to  Goncalves  et  except  for  using 
Rcooi  =  0.625  K/W,  we  estimate  Pgen  =  48.2  W  at  0.4%  efficiency.  This  discrepancy 
is  primarily  due  to  the  TEG  cold  side  temperatures — their  model  assumes  these  to 
be  325  K,  while  our  more  realistic  case  estimates  them  to  average  427  K.  Not  only 
does  this  affect  Pgen  through  Th-  Tc  (Eq.  6),  but  the  TE  properties  also  change. 

4.  Conclusion 


This  study  considers  the  temperature  load-leveling  capability  of  integrated  heat 
pipes  and  spreaders  in  a  vehicle  engine  exhaust  WHR  system  using  TEGs.  While 
simple,  this  study  highlights  important  points  to  consider  when  designing  such  a 
system: 
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1 .  TEG  properties  are  extremely  sensitive  to  temperature.  Designing  a  WHR 
system  for  a  single  operating  point  will  lead  to  very  poor  performanee  when 
applied  to  a  real  system,  as  temperature  fluetuations  drastieally  ehange  the 
TE  material  ZT. 

2.  Development  of  higher  temperature,  higher  ZT  TE  materials  will  make 
WHR  more  praetieal.  In  systems  using  eommereially  available  TE 
materials,  mueh  exhaust  enthalpy  is  wasted  beeause  the  TE  performanee  is 
so  poor  at  elevated  temperatures  (-800-1100  K). 

3.  Integrating  heat  pipes  and  spreaders  in  the  eross-flow  heat  exehanger 
system  provides  very  good  temperature  load  leveling,  whieh  holds  promise 
for  future  WHR  applieations.  The  downside  is  that  extraeting  and 
transporting  heat  from  the  exhaust  stream  to  the  TEG  hot  sides  requires  very 
low  thermal  resistanee  paths.  Inserting  multiple  heavily  finned  heat  sinks 
into  the  exhaust  stream  will  obstruet  the  flow  and  eould  lead  to  elevated 
baekpressures  and  deereased  engine  performanee. 

4.  Many  studies  in  the  literature  make  poor  assumptions  regarding  system 
operating  eonditions  and  TE  material  properties,  whieh  leads  to 
overestimated  results.  Reetifying  this  issue  will  require  investments  in  both 
improved  system-level  modeling,  aeeounting  for  the  full  heat  flow  path 
under  realistie  usage  eonditions,  and  experimental  work  to  produee 
validated  performanee  profiles  of  integrated  WHR  system  eomponents. 
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